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Abstract 

We present recent results of hadron spectroscopy and hadron -hadron interaction from the perspective of constituent quark models. 
We pay special attention to the role played by higher order Fock space components in the hadron spectra and the connection of this 
extension with the hadron-hadron interaction. The main goal of our description is to obtain a coherent understanding of the low- 
energy hadron phenomenology without enforcing any particular model, to constrain its characteristics and learn about low-energy 
realization of the theory. 
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1. Introduction 

For almost thirty years after the discovery of the y/i/r and its excitations, the so-called November revolution 
heavy hadron spectroscopy was at rest. Paraphrasing Lord Kelvin famous speech |2], by 2003 there were only two 
clouds on the horizon obscuring the beauty and clearness of the dynamical theory. In the hadron spectra these two 
clouds were on the one hand the missing resonance problem, i.e., all quark models predict a proliferation of excited 
states which have not been measured, and, on the other hand, the observation by BaBar of an open-charm meson, the 
Z)*^(2317), whose properties were quite different from those predicted by quark potential models. 

Since those peaceful days an increase on both amount and quality of experimental data has shown quite a different 
picture, far more involved and convoluted, with the observation of the well-known X(3872), the X(3940), F(3940), 
Z(3930), 7(4140), and several other states. Based on Gell-Mann conjecture ^ the hadronic experimental data were 
classified either as qq or qqq states according to S U(3) irreducible representations, nowadays this hypothesis may be 
in question. Therefore, the study of the role played by higher order Fock space components in the hadron spectra, 
allowed by the Gell-Mann classification, is an interesting issue to address. 

In this talk we give an overview of a project for getting a coherent understanding of the low-energy hadron 
phenomenology from the perspective of constituent quark models |4]. We will make emphasis on three different 
aspects. Firstly, we will address the study of meson spectroscopy in an enlarged Hilbert space considering many-quark 
components. This seems nowadays unavoidable to understand the experimental data and it builds a bridge towards 
the description in terms of hadronic degrees of freedom. Secondly, the same scheme used to describe the hadron 
spectroscopy should be valid for describing the low-energy meson-meson scattering. In other words, contributions 
to the hadron spectroscopy arising from many-quark components allowed in the Gell-Mann scheme, should also be 
recovered by means of hadron-hadron scattering using a full set of states. Finally, we will review our recent efforts to 
describe the structure of heavy baryons containing charm or bottom quarks. 

1 



/. Vijande et al. /Nuclear Physics A 00 (2012) /-[l 



2 



2. Meson spectroscopy beyond the naive quark-model 

Many-quark systems present a richer, and therefore more complicated, color structure than standard hadrons built 
with a quark-antiquark pair or three quarks. While the color wave function for ordinary mesons and baryons leads 
to a single vector, in the case of four-quark states there are different vectors driving to a singlet color state out of 
either colorless or colored quark-antiquark two-body components. Thus, when dealing with four-quark states a basic 
important question is whether we are in front of a single colorless meson-meson molecule or a compact state, defined 
as a system with two-body colored components. Whereas the first structure would be natural in the naive quark model, 
the second one would open a new area in hadron spectroscopy. 

Let us start by defining a physical channel as a four-quark state made of two quark-antiquark color singlets. In 
Ref. iHt] the formalism to evaluate the probability of the different physical channels for an arbitrary four-quark wave 
function was derived. For this purpose any hidden-color vector of a four-quark color basis, i.e., vectors with non- 
singlet internal color couplings, was expanded in terms of singlet-singlet color vectors. Such a procedure gives rise 
to a wave function expanded in terms of nonorthogonal vectors belonging to different orthogonal basis, and therefore 
the determination of the probability of physical channels becomes cumbersome. Such nonorthogonal expansions are 
common in other fields of science, like chemical physics [6], however they have not been properly discussed until now 
in the context of the multiquark phenomenology, i.e., color chemistry. We have derived the two hermitian operators 
that are well-defined projectors on the physical singlet-singlet color states | 1 l)c and | l'l'}c- Using these operators 
the probabilities for finding singlet-singlet color components for an arbitrary four-quark state I*?) are given by, 

= i :r\('i'\PQ\'i') + ('i'\QP\'i')] 

2(1 - 1,(11 I l'l'),p) L\ I ^1 / \ 1^ I /J 

^"'^[^'^'^^ ^ 2(i-i,<niPi'),p) + (1) 

where and 2(2) are the projectors over the different color vectors, | ll)c(| l'r)c)and | 88)^, ( | 8'8'),) respec- 
tively. 

The stability of a four-quark state relies on - E^q - E{Mi , M2), the energy difference between its mass and 
that of the lowest two-meson threshold. It is important to emphasize the relevance of comparing four-quark energies 
with respect to thresholds obtained using the same quark-quark interaction as in the four-quark case. Only in doing so 
one may guarantee that no spurious bound states are obtained. < indicates that all fall-apart decay channels are 
forbidden, and therefore one has a proper bound state. Ag > will indicate that the four-quark solution corresponds to 
an unbound threshold (two free mesons). As we will discuss below, between these two limits one may find molecular 
states. 

The four-body QQhh (Q stands for a heavy c or b quark and n for a light u, d or s quark) Schrodinger equation has 
been solved using two independent methods, the hyperspherical harmonic formalism ff\ and a variational approach 
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Table 1. Four-quark state properties for selected quantum numbers. The notation M1M2 If stands for mesons Mi and Mi with relative orbital 
angular momentum £. P[|33>'^(|66>'^)] and P[ | 11)^ ( I 88)^)] stand for the probability of the color components in different basis 
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Figure 1. Position probability distributions corresponding to a unbound: (S = 1, / = 1) cam, a molecular: (S = 0, / = 1) ccnn, and a bound state: 
(5 = 1,7 = 0) hbiln. 

based on generalized gaussians fs*]. All possible quantum numbers for both double charm and double bottom four- 
quark states were evaluated. Among all possible combinations, in the charm sector only one bound state with quantum 
numbers = (O)l^ has been found. We show in Table [T] a summary of the results obtained for several bound and 
unbound four-quark states in the bottom and charm sectors. One can see how, independently of the binding energy, all 
of them contain a sizable octet-octet color component. Let us first of all concentrate on the two unbound states. Using 
the formalism of Ref. |5] one can evaluate the probability of the different physical channels: Pmm, two pseudoscalar 
mesons, Pm-M', two vector mesons, and Pmm', a pseudoscalar and a vector meson. One can see how four-quark 
unbound states are represented by two isolated mesons. Let us now turn to the bound states. In contrast to the 
results obtained for unbound states, the probabilities in several of the allowed physical channels are relevant. Thus, it 
becomes clear how the bound state must be generated by an interaction that is not present in the asymptotic channel, 
sequestering probability from a single singlet-singlet color vector due to the interaction between color octets. Such 
systems are clear examples of compact four-quark states, in other words, they cannot be expressed in terms of a single 
physical channel. 

To illustrate in more detail the differences observed in the calculated four-quark wave functions we depict in Fig.[T] 
the position probability distributions defined as 

Rira, rj}) = r^rp ^ J" \Ri{ra, rp, ry)\^dry dQ.r„ dQ^^ , (2) 

where /?,(r„, r^, ry) are the four-quark radial wave functions. We present results for a unbound, a molecular and a 
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bound state, showing the position distribution for the different planes (r„, r^) = (z, x), (z,y), and (x,y). Clear differ- 
ences among them can be observed. The position distribution for the unbound case spreads in the x and y Jacobi 
coordinates up to 60 fm, while the bound and molecular systems are restricted to the region below 3 fm (molecular) 
and 1 fm (bound). In the (x,y) plane the unbound state is so widely spread that the values for the position distribution 
are three orders of magnitude lower than in the {z,y) and (z, x) cases, and therefore they will not appear in the picture 
unless artificially magnified. In the case of the molecular state a long range tail propagating in the x = y region can be 
observed contrary to the constrained values obtained for bound systems. 

In a recent investigation, the QQnh Schrodinger equation has been solved accurately using the hyperspherical 
harmonic formalism [9]. All possible quantum numbers were discussed by means of different constituent quark 
models widely used in the literature. No compact bound states were found for any set of quantum numbers nor any 
constituent quark model. Thus, independently of the quark-quark interaction and the quantum numbers considered, 
the system evolves to a well separated two-meson state. In any manner one can claim for the existence of a compact 
bound state for the ccnh system. 

Unfortunately, close to a threshold, methods based on a series expansion fail to converge since arbitrary large 
number of terms are required to determine the wave function. For this reason, we have solved the two-body Lippmann- 
Schwinger equation for negative energies using the Fredholm determinant method, looking for attractive channels 
that may contain a meson-meson molecule, that permitted us to obtain robust predictions even for zero-energy bound 
states 1 101 . To do so we have started from a physical system made of two mesons, M\ and M2, with quantum numbers 
(I)J'''^ in a relative S -wave. Then, in general, there is a coupling to the M1M2 Z)-wave and to any other two-meson 
system that can couple to the same quantum numbers (7)7^''. Thus, the Lippmann-Schwinger equation for the M-1M2 
scattering becomes 



(p,, pp- E) = Vl";f/'\p„, pp) + 2 Z P^dPyyay-M^'iPa, Py)Gy(E- Py)t';;^f/'\Py, Pp, E) , (3) 

1/ / _n ^0 



r C=o,2 



where a = fi = M1M2, y stands for any intermediate two-body state that can couple to the {T)J^'^ quantum numbers, 
f is the two-body scattering amplitude, J, I, and E are the angular momentum, isospin and energy of the two-body 
system, (aSa, (ySy, and (psp are the initial, intermediate, and final orbital angular momentum and spin, respectively, 
and py is the relative momentum of the two-body system y. To solve the scattering problem we have derived a meson- 
meson potential in the Bom-Oppenheimer approach from the basic qq interaction used to study the four-quark system 
with the hyperspherical harmonic formaUsm or the generalized gaussian variational approach. 

We have analyzed all positive-parity channels made by S -wave interacting mesons up to - 2^ in the QQnh 
sector Only the channel (/)7^ = (O)l^ is attractive enough to be bound [111 . This is the same state that was found 
to be bound in the four-quark study. As can be seen in Fig.|2]the DD* and D*D* potentials are attractive, however 
none of them is bound by itself. A bound state is obtained when one considers the coupling between the DD* and the 
D*D* systems and solves the coupled-channel problem accordingly. This mixing at short distances will reconstruct 
the octet-octet color vector component of a particular basis as shown in Table[T] 

When moving to the bottom sector one finds four more candidates for observation. They are, ordered by their 
binding energy: {I)J^ - (O)l^, (0) 0^, (1) 3" and (0) 1". The interest in these states is increasing as for the first time 
there are chances to observe such large mass exotic states in quite a near future. LHC may discover tetraguark states 
via gluon-gluon fusion due to both large number of events and their unique signature in the detectors fl2]. RHIC 
might identify hadronic molecular states by means of relativistic heavy ion collisions by employing the coalescence 



model for hadron production 111311 . If any of these suggestions gets real, the new data related to these double charm or 
bottom exotic systems will have a huge impact on heavy quark spectroscopy. 

In the QQnh case, DD and charmonium-light two-mesons systems get coupled. A few channels are slightly 
attractive (Fredholm determinant smaller than one), DD with (I)J^'^ = (0)0^^, DD with (0)1^^ and D*D with 
(0)0^^, (0)2^^, and (1)2^+. The remaining ones are either repulsive or have zero probability to contain a bound state 
or a resonance. The only bound state appears in the (7)7^'' - (0)1^^ channel as a consequence of the coupling between 
DD* and J/^'co two-meson channels. It has a binding energy in the range 0-1 MeV, and therefore in agreement with 
the experimental measurements of the X(3872). The corresponding isovector channel becomes repulsive due to the 
coupling to the lightest channel that includes a pion. Therefore, the existence of charged partners for the X(3872) can 
be discarded as well as a partner in the bottom sector iflil . 
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Figure 2. Fredholm determinant for the (/)/'' = (0) 1 ^ ccnri system. Tlie daslied line corresponds to the calculation considering only DD' singlet- 
singlet color states. The dashed-doted line stands for the results considering only the D'D' system. The solid line represents the results including 
the coupling DD' D'D' . 



Although the last analysis has been performed by means of a particular quark interacting potential, the conclusions 
derived are general. They mainly rely on using the same hamiltonian to describe tensors of different order, two- and 
four-quark components in the present case. When dealing with a complete basis, any four-quark deeply bound state 
has to be compact. Only slightly bound systems could be considered as molecular. Unbound states correspond to a 
two-meson system. A similar situation would be found in the two baryon system, the deuteron could be considered 
as a molecular-like state with a small percentage of its wave function on the AA channel, whereas the H-dibaryon 
would be a compact six-quark state. When working with central forces, the only way of getting a bound system is to 
have a strong interaction between the constituents that are far apart in the asymptotic limit (quarks or antiquarks in the 
present case). In this case the short-range interaction will capture part of the probability of a two-meson threshold 
to form a bound state. This can be reinterpreted as an infinite sum over physical states. This is why the analysis 
performed here is so important before any conclusion can be made concerning the existence of compact four-quark 
states beyond simple molecular structures. 

Lets conclude this section by emphasizing that at short distances identical quarks can recouple to different vectors 
of the Hilbert space. If any of them is not considered, a fundamental ingredient of the calculation would be neglected. 
Such effect would never happen when dealing with hadronic degrees of freedom. Therefore, the need to incorporate 
the complete Hilbert space is evident in both approaches. If the model parameters were fitted to some observables they 
would necessarily include the effect of a restricted Hilbert space and one should never use them for different quantum 
numbers without arriving to wrong conclusions, loosing thus any predictive power We can therefore conclude that 
formalisms based on four-quark and meson-meson configurations are fully compatible whenever all the relevant basis 
vectors are taken into account. 



3. Heavy baryon spectroscopy 

Heavy baryons containing a single heavy quark are particularly interesting. The light degrees of freedom (quarks 
and gluons) circle around the nearly static heavy quark. Such a system behaves as the QCD analogue of the familiar 
hydrogen bounded by the electromagnetic interaction. When the heavy quark mass niQ oo, the angular momentum 
of the light degrees of freedom is a good quantum number Thus, heavy quark baryons belong to either SU(3) 
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Table 2. Masses, in MeV, of charmed and bottom baryons. 



antisymmetric 3f or symmetric 6f representations. The spin of the light diquark is for 3f, while it is 1 for 6f. Thus, 
the spin of the ground state baryons is 1 /2 for 3f, representing the A/, and baryons, while it can be both 1 /2 or 3 /2 
for 6f, allocating 2^,, S^, HJ^, B^, Q;, and Q,*, where the star indicates spin 3/2 states. Therefore heavy hadrons form 
doublets. For example, S;, and SJJ will be degenerate in the heavy quark limit, being their mass splitting caused by the 
chromomagnetic interaction at the order l/niQ. These effects can be, for example, taken into account systematically 
in the framework of heavy quark effective field theory. The mass difference between states belonging to the 3f and 6f 
representations do also contain the dynamics of the light diquark subsystem, hard to accommodate in any heavy quark 
mass expansion. Therefore, exact solutions of the three-body problem for heavy hadrons are theoretically desirable 
because they will serve to test the reliability of approximate techniques: heavy quark mass expansions, variational 
calculations, or quark-diquark approximations. 

The combined study of Qnn and QQn systems will also provide some hints to learn about the basic dynamics 
governing the interaction between light quarks. The interaction between pairs of quarks containing a heavy quark Q 
is driven by the perturbative one-gluon exchange. For the Qnn system the mass difference between members of the 6f 
SU(3) representation comes determined only by the perturbative one-gluon exchange, while between members of the 
6f and 3f representations it presents contributions from the one-gluon and pseudoscalar exchanges. If the latter mass 
difference is attributed only to the one-gluon exchange (this would be the case of models based only on the perturbative 
one-gluon exchange), it will be strengthened as compared to models considering pseudoscalar potentials at the level 
of quarks, where a weaker one-gluon exchange will play the role. When moving to the QQn systems only one-gluon 
exchange interactions between the quarks will survive, with the strength determined in the Qnn sector, where we have 
experimental data. This will give rise to larger masses for the ground states, due to the more attractive one-gluon 
exchange potential in the Qnn sector, what requires larger constituent quark masses to reproduce the experimental 
data. 

In Ref. fl3\ the three-body Schrodinger equation was solved by the Faddeev method in momentum space with the 
constituent quark model of Ref. [ 16]. The results are shown in Table |2] compared to experiment and other theoretical 
approaches. All known experimental data are nicely described. Such an agreement and the exact method used to solve 
the three-body problem make our predictions also valuable as a guideline to experimentalists. 

It is worth noticing that the relativistic quark-diquark approximation lll9ll and the harmonic oscillator variational 
method fis"] predict a lower 3/2+ excited state for the A;, baryon. Such result can be easily understood by the 
influence of the pseudoscalar interaction between the light quarks on the Ai(l/2+) ground state. If this attraction were 
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Table 3. Ground state and excitation energies, Ail, of double charmed and bottom baryons. Masses are in MeV. 



not present for the Afc(l/2+), the Afc(3/2+) would be lower in mass as reported in Refs. (a similar effect will 

be observed in the charmed baryon spectra). Thus, the measurement and identification of the Ai(3/2+) is a relevant 
feature that will help to clarify the nature of the interaction between the light quarks in heavy baryon spectroscopy. 

Finally, we can make parameter free predictions for double charmed and bottom baryons. Our results are shown 
in Table[3] For double charmed baryons, the ground state is found to be at 3579 MeV, far below the result of Ref. 118 1 
and in perfect agreement with lattice nonrelativistic QCD li20ll . but still a little bit higher than the non-confirmed 
SELEX result, 3519 MeV |21]. It is therefore a challenge for experimentalists to confirm or to find the ground state 
of double charmed and bottom baiyons. 



4. Summary 

To summarize this talk, we have reviewed some recent results of hadron spectroscopy and hadron-hadron inter- 
action from the point of view of the constituent quark model. The constituent quark model is a powerful tool at 
our disposal to face the new issues appearing in an almost daily basis in the heavy quark sector Unfortunately, the 
downsides of such robust tools lie in the risk of misusing them. Therefore, we have tried to impress upon the reader 
the range of validity and proper formalisms of the constituent quark model. Similarly, since there are many different 
models in the literature, we have tried to present conclusions that are independent of the particularities of the model 
chosen. We have seen how the enlargement of the Hilbert space when increasing energy, that was seen to be necessary 
to describe the nucleon-nucleon phenomenology above the pion production threshold, seems to be necessary to under- 
stand current problems of hadron spectroscopy. We have also tried to emphasize that spectroscopy and interaction can, 
and must, be understood within the same scheme when dealing with quark models. Any other alternative becomes 
irrelevant from the point of view of learning about properties of QCD. 

It is expected that in the coming years better-quality data on all sectors discussed in this talk will become available 
and our results can be used to analyze these upcoming data in a model-independent way. 
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